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ABSTRACT 

We describe the initial results of a deep long-slit emission line search for 
redshifted (2.7 < z < 4.1) Lya. These observations are used to constrain 
the fluorescent Lya emission from the population of clouds whose absorption 
produces the higher-column-density component of the Lya forest in quasar 
spectra. We use the results to set an upper limit on the ultraviolet ionizing 
background. Our spectroscopic data obtained with the Keck II telescope at 
A/AAfwhm ~ 2000 reveals no candidate Lya emission over the wavelength range 
of 4500-6200 A along a 3' slit in a 5400 s integration. This null result places the 
strongest limit to date on the ambient flux of Lyman continuum photons at z ~ 3. 
Typically we attain ala surface brightness sensitivity to spectrally unresolved 
line emission in a 1 □" aperture of (0.4-0.6)xl0“^® erg s“^ cm“^ arcsec“^, and 
we search for extended emission over a wide range of spatial scales. Our 3 a 
upper bound on the mean intensity of the ionizing background at the Lyman 
limit is Juo < 2 x 10“^^ erg s“^ cm“^ Hz“^ sr“^ for 2.7 < < 3.1 (where we 

are most sensitive), assuming Lyman limit systems have typical radii of 70 kpc 
(go = 0.5, Hq = 50 km s“^ Mpc“^). This constraint is more than an order of 
magnitude more stringent than any previously published direct limit. However, 
it is still a factor of three above the ultraviolet background level expected due 
to the integrated light of known quasars at z ~ 3. This pilot study conhrms the 
conclusion of Gould & Weinberg (1996) that integrations of several hours on a 
10-m class telescope should be capable of measuring J^q at high redshift. Our 
results suggest that the integrated flux of Lyman continuum photons escaping 
from star-forming galaxies at these epochs cannot exceed twice that from known 
quasars, and that the completeness of optically-selected quasar catalogues 
must be better than 30%. We also show that it is unlikely that decaying relic 
neutrinos, if comprising the bulk of the dark matter, are responsible for the 
meta-galactic radiation held. 
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1. Introduction 

In most models of structure formation the majority of baryonic matter in the Universe 
is contained in a uniformly-distributed intergalactic medium (IGM), which is predicted 
as a product of primordial nucleosynthesis. The Gunn-Peterson constraint on the H I 
component of this IGM (Gunn & Peterson 1965; Scheuer 1965) implies that the gas must 
have been highly ionized by z ~ 5, because the flux decrement observed on the blue side of 
the Lya emission lines in the spectra of high-redshift quasi-stellar objects (QSOs) appears 
to be entirely consistent with absorption by discrete clouds along the line of sight (Steidel 
& Sargent 1987; Giallongo et al. 1994). 

The meta-galactic ultraviolet (UV) background, which is widely believed to be 
primarily the integrated light of QSOs, is thought to be responsible for maintaining the 
diffuse IGM and the Lya forest clouds in a highly-ionized state. This UV background may 
also be responsible for the ionization of metal-rich QSO absorption systems (Steidel 1990), 
the ionization of H I clouds in the Galactic halo (Ferrara & Field 1994), as well as for 
producing the sharp edges of H I disks in nearby spirals (Dove & Shull 1994; Maloney 1993; 
Bland-Hawthorn, Freeman & Quinn 1997). 

From the evolution of the QSO luminosity function, recent theoretical calculations 
for 2 ; ~ 2.5-3 indicate that the contribution of observed QSOs alone to the mean 
intensity of the UV background is Jyo(- 2 i) ~ 0.6 at the Lyman limit, i.e., at a frequency 
vq = c / 912 A (Haardt & Madau 1996, hereafter HM96). Throughout we adopt the notation 
JuQ(- 2 i) = Juo / 10“^^ erg s“^ cm“^ Hz“^ sr“^ and we assume a cosmology of go = 0.5, Aq = 0 
and hso = Hq / 50 km s“^ Mpc“^, unless otherwise stated. Indirect estimates of J^o are 
uncertain (see §^) and, given its importance, strenuous efforts have been made to directly 
measure the meta-galactic UV background. Hogan & Weymann (1987) proposed that 
long-slit spectroscopy of “blank sky” should reveal patches of fluorescent Lya emission from 
Lya forest clouds excited by the meta-galactic background. A detailed radiative-transfer 
treatment of Lya emission from thick clouds (Vhi > 10^° cm“^) photo-ionized by the 
meta-galactic radiation is given by Dinette et al. (1993). 

In a more recent paper, Gould & Weinberg (1996) show that for clouds with column 
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density iVni ^ 3 x 10^® cm“^ (optically thick to Lyman continnnm photons), ~ 50% of the 
incident energy of the ionizing backgronnd emerges in the form of Lya photons, inclnding 
those resnlting from collisional excitation. This fraction is robnst and independent of clond 
geometry, implying nniform snrface brightness. A Lya photon is absorbed and re-emitted 
within the clond nntil it is Doppler shifted by an atom with a velocity ± 4 at which 
point it can escape. A clond with a typical one-dimensional (ID) velocity dispersion of 
= 30 km s“^ wonld have a donble-peaked Lya line with a separation of 8 = 240 km s“^. 

Hence, moderate-resolntion spectroscopy {R = A/AAfwhm — 1500) of an optically-thick 
clond, along with adeqnate integration time, will yield a direct measnrement of the intensity 
of the ionizing backgronnd. Lowenthal et al. (1990) have observed the Lyman limit system 
(LLS) towards Q0731-I-653, and on the basis of a two honr exposnre spectrogram with the 
Mnltiple Mirror Telescope (MMT), place a 2 a npper limit of Jyo(- 2 i) < 200 a.t z = 2.91, on 
the assnmption that the clonds are only 1" in extent. Comparable sensitivity is obtained 
by Martinez-Gonzalez et al. (1995) throngh narrow-band imaging on the Isaac Newton 
Telescope (INT). Assnming larger clonds (~ 10"), their 3 a npper-limit is eqnivalent to 
Jv 0 (- 2 i) < 80 at z = 3.4 (see Fignre |I| and Eqnation ||). 

Moderate-dispersion long-slit spectroscopy appears the most promising techniqne with 
which to detect Lya flnorescence. From typical qnasar lines-of-sight, there are on average 
one or more optically-thick LLS per nnit redshift interval at ^ > 2 (Stengler-Larrea et al. 
1995). Hence, a single spectrogram with a slit a few arcmin in length will intercept several 
snch high-colnmn-density clonds (see §2.3|) . Imaging with narrow-band interference hlters 
(typically R ~ 100), althongh covering a larger solid angle over a smaller redshift interval, 
provides less contrast against the sky for line emission than does spectroscopy. Therefore, 
narrow-band imaging cannot achieve the reqnired sensitivity withont prohibitively long 
integrations. Similarly, the expected snrface brightness of Lya is at least an order of 
magnitnde too faint to be detected in the deepest broad-band optical imaging, i.e., the 
Hnbble Deep Field (Williams et al. 1996). Even if snfficient depth could be attained, since 
the broad-band covers a large slice of redshift space the measurement would be confusion 
limited, as the LLS clouds (expected to be ~ 10" in extent, § f^.3| ) would overlap. 

In this paper, we present the results of deep Keck H long-slit spectroscopy with the Low 
Resolution Imaging Spectrometer (LRIS, Oke et al. 1995) at moderately-high dispersion in 
the blue, 150 km s“^ full width half maximum (FWHM). In §|] we describe the long-slit 
spectroscopic observations obtained, detail the data reduction procedures and determine an 
upper-limit to J^q. We discuss the cosmological implications in §^. 
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2. Observations 

As part of another program, we have been nndertaking deep long-slit spectroscopy on 
the helds of QSOs whose spectra show damped Lya absorption systems. One field involved 
a set-np in the blue with a high-dispersion grating, and we are able to use the “blank sky” 
in the long-slit spectrogram for this investigation. The primary targets in these data are 
two actively star-forming galaxies at z = 2.81 (S2 & S3 in Warren & Mpller 1996) with 
angular separations of 12" & 21" from the Zem ~ ^^abs QSO/damped system PKS0528—250, 
which lies at the same redshift. The coordinates of the slit center were 0:2000 = 05^30™09!35, 
<^2000 = —25°03'28"2, and the position angle was set to 95?2 East of North so as to intercept 
both galaxies. The spectra of the target objects are to be presented in a forthcoming paper 
(Bunker et al. 1998). 

We obtained the long-slit spectroscopy on the night of 1998 January 20 U.T. using 
LRIS at the f/15 Cassegrain focus of the 10-m Keck II Telescope. The LRIS detector is a 
Tek 2048^ CCD with 24 /im pixels and an angular scale of 0"212 per pixel. We read out 
the CCD in two-ampliher mode. The observations were obtained using the 900 line mm“^ 
grating in Erst order blazed at 5500 A, producing a dispersion of 0.839 A per pixel. The 
reference arc lamps and sky-lines have FWHM 3-3.5 pixels, so for objects that fill the 
slit, the velocity width of a spectrally unresolved line is FWHM 150-175 km s“^. The 
grating was tilted to sample the wavelength range 4500-6200 A, corresponding to Lya in the 
redshift range 2.7 < 2 : < 4.1. The observations used the 1" by 3' long-slit. A total of 5400 s 
of on-source integration was obtained, and this was broken into three individual exposures 
of duration 1800 s that spanned an airmass range of 1.41-1.53. The telescope was dithered 
by 775 along the slit between each integration to enable more effective cosmic ray rejection 
and the elimination of bad pixels. In the R-band the seeing had a FWHM = 078-171, and 
the sky background (no moon) was 22.3 mag arcsec“^. 

Spectrophotometric standard stars G191B2B & HZ44 (Massey et al. 1988; Massey & 
Gronwall 1990) were observed at similar airmass to determine the spectral shape of the 
sensitivity function, but our observations were rendered non-photometric by intermittent 
thin cirrus. However, we were able to boot-strap our spectrophotometry to the narrow-band 
Lya line fluxes from Mpller & Warren (1993) of the two = 2.81 star-forming galaxies on 
the long slit. As we expect the Lya emitting patches to have dimensions much greater 
than our slit width and to have approximately uniform surface brightness, the photometric 
zero-points were corrected for slit losses using standard star observations with two different 
slit widths. 
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2.1. Data Reduction 

The initial stages of the data reduction followed standard procedures and were 
accomplished using IRAF|^. All the calibration and data frames were bias subtracted by 
htting a low-order function to the appropriate over-scan region. Each half of the CCD was 
then converted to photo-electrons by correcting for the gain of the corresponding read-out 
ampliher. A high signal-to-noise ratio (SNR) dark frame was then subtracted. 

The next step consisted of calibrating the pixel-to-pixel sensitivity variations of the 
CCD. Contemporaneous flat helds were obtained with a halogen lamp immediately after 
the science exposures, and these internal flats were normalized through division by the 
extracted lamp spectrum. The SNR for these flat helds was 100-300 per pixel. In order to 
remove the residual lower spatial frequency variations in the hat held due to non-uniform 
illumination of the slit, we also obtained spectral dome and twilight-sky hats. Cosmic 
rays were identihed in the hat helded spectra by iteratively calculating the noise in each 
column perpendicular to the dispersion axis and identifying discrepant pixels; those pixels 
associated with structure sharper than the point-spread function (PSF) were hagged as 
cosmic rays. 

Wavelength calibration was obtained from arc lamp spectra also taken immediately 
after the data frames. To quantify the distortions across the CCD, we traced the centroid 
of each arc-line (and the sky-lines in the data frames) at 10 pixel intervals over the length 
of the slit (800 pixels). It was found that the centroids of the lines shifted by up to 3 pixels, 
comparable to the instrumental resolution. This distortion was well ht by a quadratic 
along the columns (the spatial axis) and a cubic along the rows (the dispersion direction). 
Hence we were able to form a distortion matrix to map from a pixel space to a wavelength 
solution, leaving RMS residuals to the ht of 0.1 pixels (0.08 A). 

Rather than doing the sky-subtraction in the usual manner by htting a polynomial 
to each column (which might subtract the extended emission we are looking for), we hrst 
rectihed the sky lines using the distortion matrix. We then created a high SNR ID sky 
spectrum by averaging the rectihed rows. The slit illumination function was determined by 
collapsing the two-dimensional (2D) blank sky spectrum dispersion axis, and this function 
was used to scale the ID sky spectrum for row-by-row sky subtraction. 

The individual sky-subtracted, rectihed 2D spectra taken at diherent dither positions 


^The Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 
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along the slit were then registered and averaged, with a percentile clipping algorithm used 
to reject cosmic rays that had escaped previous cuts. A mask was used to flag the known 
bad pixels, which were excluded from the averaging. The individual 1800 s integrations 
were sufficiently long for our data to be background (rather than read-out noise) limited 
at all wavelengths. The measured noise agreed well with the Poissonian statistics of the 
sky background, demonstrating that any additional noise introduced by our data reduction 
techniques was insignificant. 


2.2. Analysis 

After rectification and sky-subtraction, we searched for extended Lya line emission by 
applying various filters along the spatial axis of the combined 2D spectrum. We convolved 
the 2D spectrum with a highly-elliptical 2D Gaussian oriented parallel to the slit. For 
the minor axis (along the dispersion direction), the width of the kernel was set to be 
comparable to the spectral resolution. We varied the length of the major axis to search for 
signal on any characteristic scale length between 1" (our seeing disk) and 180" (the length 
of the slit). For the anticipated range in Wq; fhe signal per resolution element is not large, 
so this smoothing is essential. We repeated the smoothing with simple boxcar filtering of 
various lengths along the slit axis, followed by convolution with a Gaussian of comparable 
extent to the seeing disk. This second approach is optimal if the Lya clouds do indeed 
have uniform surface brightness (Gould & Weinberg 1996). For both techniques, visual 
inspection revealed no significant extended single-line emission. 

The only Lya line-emission sources detected were the target galaxies, PKS0528—250 
SI & S2. These have also been observed in continuum imaging (Mqller & Warren 1998), 
with rest-frame equivalent widths of Wq = 88 A & 46 A respectively. This implies that Lya 
arises from an internal UV source (star formation or AGN activity). Gharlot & Fall (1993) 
calculate Wq = 80-200 A for a star forming galaxy, depending on the initial mass function 
(IMF). However, selective extinction of this line through resonant scattering and suppression 
by dust can greatly reduce the equivalent width of Lya in star forming galaxies (e.g., Steidel 
et al. 1996). In stark contrast, Lya produced by external photo-ionization should have a 
much larger equivalent width, as the continuum arises mainly from two-photon decays. 
Hence, the rest-UV continuum around Lya provides excellent discrimination between local 
OB stars and the diffuse meta-galactic background as the source of the Lyman continuum 
photons. Any source with a detectable continuum can only yield an upper-limit on as 
it will likely be dominated by local effects. 

Our long slit also intercepted a [O II] A 3727 A emission line aX z = 0.424, 6175 
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East of S3, associated with a V = 22.8 mag galaxy. We detect modest continuum 
(Wo[0 II] = 100 A ± 10 A), and the redshift is conhrmed with [Ne III] 3869 A. The [O II] 
doublet is well resolved, so our resolution is better than 221 km s“^. 

Our search techniques failed to reveal the signature of Lya fluorescence from 
optically-thick clouds. A lack of Lya emission constrains the incident UV flux. In order 
to test the validity of our search methods and to place meaningful upper limits on Ji^q, we 
examined the recoverability of artihcial emission lines added to our data frames. These 
synthetic lines simulated clouds of various sizes, and were convolved with the seeing disk 
and the instrumental spectral resolution. For those with shorter spatial extent, less than 
5", we were able to recover 90% of the simulated spectrally unresolved lines with total 
SNR = 3. However, for larger artificial clouds, the recoverability was slightly worse for the 
same integrated SNR (but lower surface brightness). This is because slit illumination effects 
become more dominant on these larger scale lengths, and the spectra of the target objects 
introduce confusion. For SNR = 3, the recoverability falls below 60% on scales exceeding 1'. 

Our typical noise in a 1 □" aperture corresponds to a Lya surface 
brightness of /Lya(-i8) = 0.4-0.6 (1 a). We adopt the convention /Lya(-i8) = 

-^Lya / 10“^® erg s“^ cm“^ arcsec“^. The quoted noise is for ‘clean’ regions of the 
night sky, devoid of strong lines, and the sensitivity declines slightly at the shorter 
wavelengths of our spectral range (4500 A < A < 6200 A). These limits are for a spectral 
extraction width of 4 pixels (220 km s“^). This is optimal for a ID velocity dispersion of 
= 30 km s“^, comparable to that measured for Lya forest clouds at z ^ 2.3 (Kim et al. 
1997). The line widths may in fact be greater due to coherent velocity structure (shears 
and bulk flows). For > 20 km s“^, the twin peaks of the Doppler-broadened resonant 
prohle become resolved, and the SNR falls proportional to We have corrected our 

limits for Ab = 0.05 mag of foreground Galactic extinction (Burstein & Heiles 1984). The 
recent COBE/IRAS dust maps of Schlegel, Finkbeiner & Davis (1998) indicate a greater 
extinction {Ab = 0.2 mag). 

Our 3 a upper limits on the Lya surface brightness as a function of ID velocity 
dispersion and angular cloud size, 6*, projected along our l"-wide slit are: 

/ Q \-V2 ^ 1/2 

/Lya(-18) = (0.4 - 0.6) f — j V 30kn^s~V ^ 


for 1" < 6 < 60", 2.7 < z < iA and ^ 30 km s 



2.3. Upper Limits on Uq 


We use our results to place an upper limit on Uq as a function of redshift and the size 
and velocity dispersion of the clouds. As this is derived from the limiting surface brightness 
attained, the constraint on the meta-galactic Lyman continuum flux is independent 
of Hq, qo and Aq, and is only subject to the (1 + z)‘^ cosmological surface brightness 
dimming. The measurement of W from Lya fluorescence is immune to the effects of in 
situ dust; the extra-galactic Lyman continuum photons only penetrate clouds to a depth 
of Nm ^ 10^® cm“^ corresponding to Ay ~ 5 x 10“^ mag for a Galactic gas-to-dust ratio. 
Gould & Weinberg (1996) calculate that only 4% of Lya photons are absorbed as they 
random-walk their way out of the nebula from this depth. 

An observed Lya surface brightness /Lya(-i8) would correspond to a meta-galactic 
background of: 

Vo(-2i) = 4,76 (^) (1^) (2) 

where r]E is the fraction of energy absorbed from the ionizing background that is emitted 
in the form of Lya photons (Gould & Weinberg 1996). For clouds with column depths 
exceeding Nm 3 x 10^® cm“^, this fraction is robust (p^ ~ 0.5), but falls to rjE ~ 0.4 
for Nm ~ 10^® cm“^. Note that our observations constrain directly the flux of Lyman 
continuum photons and not Wo (the energy density at the Lyman edge). Gonversion to Wo 
assumes a particular form of the spectrum. Specihcally we adopt W = Wo (^/^o)~“ so that 
W / (hpiu) dv = Wo / (hpicr), where hp\ is Planck’s constant. The HM96 spectrum has 
a slope a = 1.73. 


The greatest uncertainty on our upper limits is the size of the clouds. From a study 
of the foreground galaxies responsible for Mg II AA 2796, 2803 A QSO absorption systems 
at 0.2 < z < 1.6, Steidel & Dickinson (1995) estimate the linear extent of LLS to be 
’"LLS ~ 70 kpc. For a randomly-distributed population of spherical clouds of radius 
tlls where the angular size of the clouds greatly exceeds the slit width (1" ~ 7 h'^Q kpc at 
2 ; ~ 3), the expectation value for projected length of intersection D along our spectroscopic 
slit is: 


D 


TrrLLS ^ /8 

2 1^3 



rLLS = (110 ± 14) 


/ ^LLS \ 

V70 kpc/ 


hj kpc 


(3) 


This corresponds to an average angular projection of 6 = 15" ± 2" for go = 0.5 and 
9 = 10" ± 1" for go = 0.1 at 2 ; ~ 3, assuming no evolution in the typical size of LLS clouds. 
In calculating our limits on Wo, we conservatively adopt an average projection of 10" along 
the slit, consistent with the smallest likely size. 


The limits we place on Wo as a function of redshift are depicted in Figure |l|. Over 
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the range 2.7 < z < 3.1, a wavelength region for Lya where the sky spectrnm is relatively 
featnreless, onr sensitivity to J^o is approximately flat. This is dne to a combination of 
LRIS response, which increases towards the red, and Lya line energy flnx which decreases 
as {1 + z)~^ throngh cosmological snrface brightness dimming. In fact, onr grating 
set-np (dictated by the demands of another project) was slightly red-ward of what may 
be most desirable for low-airmass and dark-time observations. The efficiency cnrve for 
LRIS-l-900 line mm“^ grating wonld snggest that ~ 2.5 is the optimal hnnting gronnd 
(see Fignre |l|). 

Our 3 cr upper limit on the meta-galactic flux over the range 2.7 < z < 3.1 where we 
are most sensitive is given by: 

U-'n) = 2.0 (A) (^51^) (w) ^ “PP" 2'P < " < 3-1) 

(4) 

where 0 is the angular extent of the cloud, and we have assumed a resonantly-broadened line 
prohle from a cloud with velocity dispersion ay. If the optically-thick clouds are as small as 
1", the 3 a upper limit is only Jiyo(-2i) < 6.3, which is still a signihcant improvement over 
the previous limit of Ji/o(- 2 i) < 200 at 2 a from Lowenthal et al. (1990). At higher redshifts, 
our upper limits on J^q are: 


= 3.0 (A j (^) (w) " “PP” 

(5) 

= 4.0 ((^) (w) PPP“ 

( 6 ) 

A concern is that we do not know a priori that our slit actually intercepts any 
optically-thick clouds. However, the incidence in QSO spectra of such clouds with column 
densities exceeding Ahi is: 


dn 

dz 


0.9 


101V 


(^) 


(7) 


per unit redshift (HM96). Table 1 gives the likelihood that one or more optically-thick 
clouds intercept our slit, where we assume a Poissonian geometric distribution of clouds with 
a typical projection of 10". Over the redshift interval 2.7 < 2 ; < 3.1 where our constraints 
on J^Q are greatest, there should be on average ~ 1-2 such systems per arcmin surveyed. 
Therefore, we can reject the possibility that no 2.7 < 2 ; < 3.1 optically-thick absorbers lie 
on our 3' slit at the > 97% conhdence level. 
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3. Discussion 

Given the null results of our search for line emission, we consider the cosmological 
implications of our new upper limits on the ambient Lyman continuum background at high 
redshift. 


3.1. QSOs and the Proximity Effect 

The measured decrease in the number of Lyo-absorbing clouds induced by the UV 
radiation held in the neighborhood of a QSO (the proximity effect) provides an estimate 
of Jyo at high redshift, independent of Hq. A recent determination for 1.7 < 2 ; < 3.8 gives 
a mean intensity of Jvo(- 2 i) = 0.3-0.8 (Espey 1993; Lu et al. 1991). These limits are a 
factor of 2.5-7 below our detection threshold. However, the proximity effect is subject to 
several systematic uncertainties; J^q could be underestimated if the luminosity of quasars is 
highly variable, or if most of the quasars with good spectra (naturally the brightest ones) 
are magnihed by gravitational lensing, or both. From our results, the global boosting of 
the quasar luminosity function cannot exceed a factor of three. In addition, the number of 
clouds near quasars might be enhanced because of clustering. If this is not accounted for, it 
will lead to an underestimate of the ionizing photon flux. 

The QSO luminosity function and its evolution, as derived from multi-color optical 
surveys, have been used to estimate the integrated QSO contribution to the UV background. 
A hrm lower-limit over 2 < 2 : < 3 is > 0-23 (Rauch et al. 1997). Accounting for the 

opacity of H and He associated with intervening Lya clouds, Mieksin & Madau (1993) and 
HM96 estimate that the mean intensity of the diffuse radiation held from QSOs increases 
from ~ current epoch to 0.6 at z = 2.6-3.0. This is a factor of three 

below our current upper limits at high redshift. 

On one hand, this estimate of is unknown at the 50 per cent level because of 
uncertainties in the evolution of the Lyo; forest clouds (Meiksin & Madau 1993). On the 
other, this is an underestimate because some quasars at these redshifts may be missed 
through foreground dust extinction in high-column-density absorbers (Fall & Pei 1993). 
Our upper limit of Juo{- 2 i) < 2.0 indicates that at the very most, two out of three QSOs at 
z ze t!) are missed through absorption by foreground objects. This is at the level where we 
can begin to constrain the models of Fall & Pei (1993), who predict that between 10% and 
70% of 2 : ~ 3 QSOs are absent from optical samples. 
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3.2. Constraints on High- 2 : Star-Forming Galaxies as the Source of Lyman 

Continuum Photons 

It has been suggested that the young, hot OB stars in actively star-forming galaxies 
at high redshift may contribute signihcantly to the meta-galactic UV background (e.g., 
Songaila, Cowie & Lilly 1990). It is likely that much of the Lyman continuum radiation is 
absorbed by H I in the parent galaxy, i.e., the H II regions are ionization and not density 
bounded. Our upper limit on at z ^ 3 implies that the number of escaping Lyman 
continuum photons from star-forming galaxies can be a factor of no more than two greater 
than that contributed by quasars, assuming ~ (HM96). 

If we take the global star formation density from known high-redshift galaxies 
(e.g., Madau et al. 1996; Madau 1997), then we can constrain the escape fraction 
(/esc) of UV photons at wavelengths short-ward of the Lyman limit for these 
galaxies. From the HDF, the co-moving luminosity density at A^est ~ 1500 A is 
p(1500 A) = 2.1 X 10^® hso erg s“^ Hz“^ Mpc“^ at x ~ 3 (Dickinson 1998). For 
a Salpeter IMF (with 0.1 Mq < M* < 125 Mq), this non-ionizing UV luminosity 
translates to an average star formation rate (SFR) per unit co-moving volume of 
Pg^( 2 : ~ 3) = 0.026 hso Mq yr“^ Mpc“^, assuming no obscuration. However, from recent 
near-infrared observations of nebular emission lines, Pettini et al. (1997, 1998) estimate 
that the UV continuum at x ~ 3 is extinguished by H 1500 ~ 1^2 mag, implying that the 
total star formation rates are a factor of 3-6 higher. This is broadly consistent with the 
recent HDF sub-mm detections by SCUBA/JCMT (Hughes et al. 1998) corresponding to 
PsFR = 0.1 hso Mq yr“^ Mpc“^ at x ~ 3. Applying the Calzetti (1997) extinction curve 
derived from local star-burst galaxies to the x ~ 3 population yields a similar rest-UV 
extinction (a factor of 2-6; Calzetti 1998). 

A constant star formation rate of 1 Mq yr“^ will produce 1.44 x 10^^ s“^ Lyman 
continuum photons, from the most recent version of the spectral evolutionary code of 
Bruzual & Chariot (1993) with the same Salpeter IMF. Hence, the x ~ 3 galaxies produce 
Lyman continuum photons at a volume-averaged rate of: 

/ 10O.4yii5oo \ /X-i-;2\3 

hLyC = 2.45 X 10“^° (- - -j f—-—j hso photons s“^ cm“^ (8) 

where fiLyC is in physical rather than co-moving co-ordinates. The mean free path of 
a Lyman continuum photon, Al, is inferred from Madau (1992), who shows that unit 
optical depth is achieved for Ax 0.2 at x = 3, corresponding to a physical length of 
Al ^ 40 h'^Q Mpc (i.e., the ionizing flux is largely local, and evolutionary effects can be 
neglected as Lyman-continuum photons from galaxies at higher redshifts are either severely 
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absorbed or redshifted to energies below the Lyman edge). Therefore, the contribution of 
star-forming galaxies to is given by: 


tSFR 

"'i/O 


/e: 


hp\Ci 

dTT 


^LyC 


(9) 


for an ionizing spectrum of the form J^, = Jvo{v/vq) which can be expressed as: 


tSFR 

■Jv0{-21) 




Al \ 
40 h^Q Mpc ) 


( PSFR \ 

V0.026 ho) 


'X00-4^15Oo\ z 


( 10 ) 

Therefore, the dust-corrected SFR per unit co-moving volume at ~ 3 is equivalent to 
Juo(- 2 i) = (3 — 6) /esc, if the rest-UV is extinguished by factors of 6-3. Hence, our 3 a 
limit of Juo{- 2 i) = 2.0 constrains /esc < 25%-50%, assuming QSOs contribute 
(HM96). However, Meurer et al. (1997) argue that in situ dust obscuration leads to a 
much greater underestimate of the integrated star formation rate, and that the rest-UV 
extinction correction should be a factor of ~ 15 rather than 3-6. If this is so, our constraint 
is signihcantly tighter (/esc < 10% at 2 : ~ 3), comparable to but still somewhat above 
the best limits on the Lyman continuum escape fraction from galaxies today (/esc < 1%, 
Deharveng et al. 1997). 


3.3. Gunn-Peterson Constraint on the Diffnse IGM 

The absence of an absorption trough on the blue side of the Lyo emission line in 
the spectra of high-redshift QSOs represents the best limit on the density of a neutral, 
smoothly-distributed IGM (Gunn & Peterson 1965; Scheuer 1965). A measurement of 
the ionizing meta-galactic flux, when combined with limits on the baryon density from 
primordial nucleosynthesis and the Gunn-Peterson optical depth (Jenkins & Ostriker 1991; 
Giallongo, Gristiani & Trevese 1992) measures what fraction of gas has collapsed into 
galaxies and discrete clouds. Our upper limit on J,^o requires that baryonic matter in 
the diffuse IGM comprises Oigm < 0.016 Ocrit of the closure density at 2 : ~ 3 (from 
Equation 2 of Songaila, Gowie & Lilly 1990). We assume the mean optical depth of the IGM 
is Tz ~3 < 0.02 (Steidel & Sargent 1987). From standard Big Bang nucleosynthesis (Walker 
et al. 1991), our results imply that the fraction of baryons residing in the diffuse IGM at 
2 : ^ 3 is at most 60% / hioo (where hioo = Hq / 100 km s“^ Mpc“^), i.e., no constraint for 
low values of Hg. We are undertaking deeper integrations which should improve this upper 
limit by a factor of three. 
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3.4. Decaying Relativistic Neutrinos and the Meta-Galactic UV Background 

The recent Super-Kamiokande detection of neutrino oscillation (Fukuda et al. 1998), 
although providing only a mass difference between neutrino flavors rather than an absolute 
measurement of mass, has rekindled interest in Hot Dark Matter (HDM). We can use 
our upper limits on J^o to constrain any relic decaying neutrino flux (e.g., Sciama 1990). 
According to this theory, relativistic neutrinos make up the bulk of the critical matter 
density (Dy ~ 0.9Dcrit), and their decay photons can ionize hydrogen (i.e., > 13.6 eV) 

but not helium. Adopting the values from Sciama (1995, 1998) for the decay photon energy 
[E-y ^ 13.8 eV), the neutrino lifetime ~ 2 x 10^^ s) and number density at the current 
epoch (W = 112.6 cm“^) would result in a UV flux at z ~ 3 equivalent to Jjyo(- 2 i) = 4.5 
(Figure 2). This would correspond to a 6 a signal in our spectrum over 2.7 < z < 3.3. Our 
failure to detect such Lya emission rules out with high conhdence a Dm = 1 cosmology in 
which the relativistic decaying neutrinos of Sciama (1998) form the bulk of dark matter. 
Furthermore, in a Mixed Dark Matter (MDM) scenario, an HDM component in the form of 
F'.y = 13.8 eV Sciama neutrinos can comprise at most VLy < 0.4Dcrit- 

Our observations constrain the combination of parameters {E^ — 13.6 eV) / Ho), the 
excess photon energy above the Lyman limit available for photo-ionization. If the neutrinos 
do account for most of the critical matter density in a flat Universe (i.e., VLy ^ 0.9 Dcrit, 

Dbaryon ^ O-l D^rit)) tlieU 

{E^ — 13.6 eV) = 0.08 x -— eV (3 a upper limit) (11) 

2 X 10"''^ s 

Therefore, the energy of the decay photons can be at most only 0.6 per cent above the 
ionization potential of H I, which would require much hne-tuning. 


4. Conclusions & Future Work 

Our moderate-dispersion long-slit spectroscopy with LRIS/Keck has placed a 3 a upper 
limit on the ambient UV background at ~ 3 equivalent to a flux at the Lyman limit of 
Jvo{- 2 i) < 2.0, assuming optically-thick clouds with dimensions ~ 10". We anticipate even 
greater sensitivity in the blue for observations at lower air mass and unaffected by cirrus. 
Although this direct limit is almost two orders of magnitude lower than any previously 
published, it is still a factor of three above the expected contribution of known QSOs at these 
epochs. However, it does enable ns to say that the completeness of optical QSO catalogs 
is better than 30%, and that the contribution to Juq{z ~ 3) from star-forming galaxies 
cannot exceed twice that from known quasars. We constrain the escape fraction of Lyman 
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continuum photons from star-forming galaxies at these redshifts to be /esc < 10%-50% 
(assuming extinction corrections of 15-3 at Arest ~ 1500 A). Although our relatively short 
(5400 s) integration in this pilot study is insufficient to reach the lower limit on the quasar 
contribution to we are able to exclude the possibility that relic neutrinos comprise the 
bulk of dark matter, as we would detect Lya fluorescence from the flux of ionizing decay 
photons. 

The size of optically-thick QSO absorption systems is largely unconstrained (§ p.3|) , and 
their morphology is entirely unknown because quasar sight-lines yield only one-dimensional 
information. Some work has been done using two lines-of-sight provided, for example, by 
a pair of lensed quasar images (e.g., Dinshaw et al. 1994). The single example of multiple 
sight-lines to a damped system indicates a size greater than 16 h. 50 ^ kpc (Briggs et al. 1989). 
Without strong constraints on the sizes of the high-column-density systems, it is impossible 
to differentiate between the scenario where outer galactic halos are directly responsible for 
the absorption or whether the culprit is more extended diffuse gas. This may be associated 
with a cluster, possibly with a fragmented, hlamentary topology (as suggested by recent 
SPH simulations, e.g., Katz et al. 1996). A detection of fluorescent Lya would reveal the 
size and morphology of the optically-thick clouds, which is a primary clue to understanding 
their origin and relation to the galaxies which are forming at these epochs. 

The phenomenal light-grasp of Keck plus the new insights afforded by Gould & 
Weinberg (1996) make it possible to design an experiment with an extremely high likelihood 
of detecting fluorescent Lya emission from high-column Lyman limit absorbers. This will 
enable the hrst direct measurement of the size of these high-redshift H I clouds, and the 
determination of Jj,o at 2 ; ~ 2.5-3. This initial study has shown that, for a dark-time 
integration of ~ 20 hours in good conditions, the sensitivity of LRIS/Keck should be 
sufficient to detect the Lyman continuum background at > 5 a, even if the meta-galactic 
flux has no major contribution other than the known high-redshift quasars. We are 
embarking on an extensive program of LRIS/Keck observations to obtain long-slit spectra 
of several hours. By using a Liter of width AA/A ~ 10%, we can obtain several parallel 
long-slits on the same spectrogram, thereby covering a large solid angle while concentrating 
on those redshifts {z ~ 2.5) where our sensitivity to Ji^ is greatest (Figure |^). For the hrst 
time we will set strict quantitative limits on the inventory of sources of UV radiation at 
high redshift. 


We are indebted to Hy Spinrad and Steve Warren for allowing us to use their 
spectroscopic data for this blank-sky search. Our thanks to Daniel Stern, Arjun Dey 
and Adam Stanford for their help during the observing run. The software for rejecting 
cosmic ray strikes was written by Mark Dickinson. Douglas Finkbeiner kindly provided us 
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Table 1. 


Upper Limits and Detection Statistics 



> 10^® 

Nm / cm 2 

> 3 X 10^® 

Ve 

0.4 

0.5 

Ju 0 {- 21 ) 3 (T limit 

2.5 

2.0 

of clouds / arcmin 

2.0 

1.2 

P('^cloud ^1) 

99.8% 

97.3% 


Note. — For the redshift interval 2.7 < < 3.1 where our limits on J^q are 

most stringent. 
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Fig. 1.— Our 3 a upper-limits on the flux density at the Lyman edge, Juo{- 2 i) = 
/ 10“^^ erg s“^ cm“^ Hz“^ sr“^, as a function of redshift (solid line). We assume a 
meta-galactic UV spectrum of the form J(z/) oc with a = 1.73 (HM96). The right axis 
shows the flux of Lyman continuum photons (/“ Jv/hp\u du), which is measured independent 
of a. The region above the heavy solid curve is excluded at the 90% confidence level. Spikes 
of reduced sensitivity occur at redshifts where Lya coincides with the sky lines in our 1.5 hour 
spectrogram. These limits assume that a long 1" slit covers ~ 10 □" of an optically-thick 
cloud (iVni 3 X 10^® cm“^) and that the cloud velocity dispersion is ay ^ 30 km s“^. Also 
shown are the previous 3 a upper-limits from the INT (Martinez-Gonzalez 1995) and the 
MMT (Lowenthal et al. 1990), converted to the same assumed projected cloud size. We 
plot our anticipated sensitivity for a 20 hour exposure with LRIS-l-900 line mm“^ at low 
airmass in dark time (dotted curved). The shaded region is the lower-limit on from 
the proximity effect (Espey 1993). The estimated contribution to from the luminosity 
function of known high-;^ QSOs is also shown as the dashed line (HM96). The equivalent 
Lyman continuum flux from decaying relic neutrinos is plotted as a dot-dash line, adopting 
the parameters of Sciama (1998) and assuming these neutrinos form the bulk of the dark 
matter in an Gm = 1 Universe. Such a scenario is strongly ruled out by the null results of 
our survey. 
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Fig. 2.— The number of sources intercepted above a given SNR threshold (SNR > 5) in a 
10-20 hour spectrum as a function of the number of slits. We perform this calculation for 
the LRIS spectrograph on Keck with the 900 line mm“^ grating set to Acent = 4047 A. Given 
the redshift distribution of Lya absorbers and our limits on Jj,o as a function of wavelength 
(Figure Q), we are most sensitive to clouds over a relatively limited range of redshifts. This 
means that pixels in the detector array are better used in the spatial dimension rather than 
the spectral dimension. Having multiple slits trades wavelength range for area covered. 
Since LRIS is equipped with a slit mask mechanism we can insert several parallel long slits 
of length 7' at the focal plane. A custom narrow band hlter (AA/A ~ 10%) would be 
used to prevent spectral overlap. This hgure shows the increase in the efficiency by using 
this multi-slit Lya search technique. We conservatively assume that the number of clouds 
per unit redshift interval scales as (1 +z), which is correct for damped Lyman alpha systems 
(Lanzetta, Wolfe Sz Turnshek 1995). For lower-column systems, dn/dz oc (1 + z)'^, where 
7 ~ 1.5 (Stengler-Larrea et al. 1995). Adoption of this steeper distribution would favor a 
multi-slit search even more. 




